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ABSTRACT 


An experiment was conducted which measured the Radio Fre¬ 
quency (RF) radiation from the PHERMEX accelerator, capable of 
30 MeV and 600 A. This was accomplished by placing TEM horn 
antennae at varying angles from the path of the electron beam. The 
signals received by the antennae were then recorded by using a Digi¬ 
tizing Camera System (DCS). Measurements were taken of the radia¬ 
tion from propagating and non-propagating beams, beams with energy 
above and below Cherenkov threshold, and beams with varied cur¬ 
rents. The captured RF signals and their corresponding frequency 
spectra were then analyzed. This analysis showed that the radio fre¬ 
quency radiation from the beams below the Cherenkov threshold 
contained primarily transition radiation; when above, diffracted 
Cherenkov radiation was observed. Non-propagating beams produced 
larger-angle radiation and had less definition in their spectrum. All 
electric fields measured were proportional to the beam current. 
Lastly, the electron beam pulse width amd separation were determined 
by both the received signals and their spectrum. 
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I. INTRODUCTION 

A. BACKGROUND 

1. Historical Look at Cherenkov Radiation 

The following history of Cherenkov radiation discovery and 
development is based on Jelley’s introduction (Ref. 1] and follows the 
summary by O’Grady [Ref. 2 ]. 

As early as 1910, the effects of Cherenkov radiation were 
observed by Mme. Curie. Her bottles of concentrated radium solution 
emitted a pale blue light. This phenomenon was subsequently partially 
Investigated by Mallet during the years between 1926 and 1929. How¬ 
ever. it was not until 1938 that Cherenkov commenced an exhaustive 
series of experiments to fully understand this phenomenon. During this 
time, a theory was proposed by Frank and Tamm which had excellent 
agreement with Cherenkov’s experimental results. Since their work, there 
have been numerous theories and experiments concerning Cherenkov 
radiation detection and Cherenkov radiation applications. 

By Frank and Tamm’s definition (Ref. 3], Cherenkov radiation 
takes place when a charged particle moves at a constant velocity which is 
larger than the phase velocity of light in a transparent medium. Electrons 
can be accelerated to speeds greater than the velocity of light in air and 
thus produce Cherenkov radiation. Frank and Tamm’s theory, however, 
contains a number of simplifying assumptions. Two of the more impor¬ 
tant assumptions are that the medium is an unbounded continuum and 
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the track length is infinite [Ref. 1]. Applying these assumptions to most 
accelerator applications is unrealistic. When an electron bunch leaves 
the vacuum of an accelerator, another medium (such as air) is encoun¬ 
tered. The electron beam steadily loses energy due to ionization of air 
molecules. Therefore, the beam has a finite length and propagates in a 
bounded medium. 

When the two previously made assumptions are violated, two 
similar types of radiation occur—diffracted Cherenkov radiation and 
transition radiation [Ref. 4]. Diffracted Cherenkov radiation is caused by 
the finite beam interaction length and transition radiation is caused by a 
boundary where tlie index of refraction of the medium changes. Another 
type of radiation closely related to transition radiation is diffracted 
transition radiation [Ref. 5|. In this case, radiation is produced when a 
charged particle passes through an aperture or travels near any interface 
between differing media. The circular beam port of an accelerator is an 
obvious source of diffracted transition radiation. However, since the 
circular beam port dimensions are constant, the resulting radiation will 
merely be added to the transition radiation and considered the electro¬ 
magnetic pulse (EMP). 

Since the major effects of this experiment concern both EMP 
and diffracted Cherenkov, a more-detailed summary for each is contained 
in the following two sections. The last section in this chapter addresses 
the similarities and differences of both. 
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2. Cherenkov and Dlffiracted Cherenkov Radiation Theory 

As previously stated, Cherenkov radiation is produced whenever 
a charged particle travels through a dielectric medium faster than the 
phase velocity of light in the same medium. The underlying reason radia¬ 
tion forms is that the charged particle temporarily excites the medium. 
The ionized molecules briefly radiate as a dipole. Of particulai' interest is 
the Cherenkov angle Oc- When the interaction length is Infinite, the radia¬ 
tion propagates only at this angle. The Cherenkov angle can be calcu¬ 
lated using 

cosGc*^ (1) 

where c is the speed of light in a vacuum, n is the refractive index of the 
medium, and v is the speed of the particle. A more-detailed description of 
Cherenkov radiation and a thorough derivative of the fundamental equa¬ 
tion for output of radiation Is given by Jelley [Ref. 1]. However, since 0c 
occurs at approximately one degree, the diffracted Cherenkov radiation 
will be more applicable to this experiment. 

Measurements by Neighbours, Busklrk, and Saglam [Ref. 6] in 1983 
showed that the Cherenkov cone angle Is broader as the Interaction 
length is decreased from infinity. This, of course, is diffracted Cherenkov 
radiation. Neighbours. Busklrk, and Saglam [Ref. 7] derive an equation 
for the total coherent power per unit solid angle, radiated at the fre¬ 
quency V by a periodic charged particle beam traveling a finite distance L 
at constant veloclly. This equation is 
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w(v.k)*v2QR® 


(2) 


where Vo is the fundamental frequency of the beam generator and where 
the constant Q is 



where \x is the permeability of the medium and q is the change in the 
electron bunch. The radiation function R is given by: 

R = W7tTisin9I(u)F(k) (3) 

where: 

e a angle between emitted radiation and beam path 

ri = L/X 

L = beam interaction length 

X = radiation wavelength 

F(k) = dimensionless form factor 

I(u) = sin u/u = diffraction function 
where u is: 


u = Taifcose^ - cose] (4) 

A graph of equation 3 with the appropriate constants, for this experiment 
yields Figure 1. 'rhis envelope function determines the relative output 
power for a given angle from tlie beam path. 


4 








Figure 1. Envelope Function for Beam Parameters 

Two additional comments concerning diffracted Cherenkov radi¬ 
ation are: (1) “the power emitted for coherent radiation Is proportional to 
the source of the beam current” [Ref. 81. (2) Neighbours, Buskirk, and 
Maruyama [Ref. 9] also predict that “the oscillations of the frequency 
distribution of the radiation are of Increasing “frequency” as the angle of 
propagation moves toward the back direction.” To better understand this 
statement, one should view Figure 2 (obtained from Reference 9). From 
Figure 2, the predicted pattern is cleai'. 


5 








:.so -r 


2.00 -3 
O : 

Z - 

UJ 


C 

UJ 

02 


«.• Z. 


1.C0 -] 


> : 

5 0.5? •: 

UJ : 

: 


/* 


Ootnta Line . / 
nJ - 0.9 . ./ 

// 

// 

/ / Solid Line 
, / ni ■ 






I I t I I l■^■^ r 'l M I I 1 i‘ I I 1 •^"r r 

20 60 90 120 150 ISO 

ANGLE (deg) 


Figure 2. Frequency Shifting with Angle 

3. Transition Radiation Theory 

While Cherenkov radiation occurs when the velocity of the par¬ 
ticle is greater than tlie speed of light in that medium, transition radia¬ 
tion can occur whether the beam energy is above or below the Cherenkov 
threshold energy. "When the electron passes through a boundary between 
two different media transition radiation is produced. Frank and Ginsburg 
[Ref. 10] explain the nature of this phenomenon by considering the 
transition of an electron from a vacuum into a perfect conductor. While 
the electron moves in the vacuum, there is no radiation. However, the 
field in the vacuum is equal to the field of the electron and its image 
moving toward it. Once the electron exits the vacuum, both the electron 
and its image cease to exist from the vacuum field perspective. Therefore, 
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radiation must be given off in much the same way as it is emitted when 
an electron suddenly stops. 

Work done by Buskirk and Neighbours [Ref. 4] develops a useful 
equation by which the electric and magnetic field strengths can be evalu¬ 
ated. The maximum magnetic field is given by 

B = i - 

^ R “1-pcose 

where 0 is the angle between the beam path and detection point, R is the 
distance between the field point and the beam port, and lo is the current 
of the accelerator. The constant P is defined to be 



where c is the velocity of radiation in the median, Co is the velocity of 
light in free space, and n is the index of refraction of the medium. For 
this experiment, the transition will be from vacuum to air as the electron 
bunch leaves the accelerator. The electric field can therefore be written as 
[Ref. Ill 


E = SS*. _ I _ (6) 

4k R “1-Pcos0 

where cdq Is the penneabillty of free space. Figure 3 [Ref. Ill shows the 
relation between transition energy emitted per angle as a function of 
beam energy. 
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Figure 3. Transition Energy Emitted per Unit Angle 

Diffraction radiation during this experiment will be measured 
and treated as additional transition radiation. 

4. Comparison of Diffracted Cherenkov and Transition 

Radiation 

The similarities between Cherenkov and transition radiation 
were summarized by Busklrk, Mack, and Neighbours [Ref. 12]. First, 
both radiation fields are polarized, with the electric field in the plane of 
the observer and the beam path, and the magnetic field perpendicular to 
that plane. Also, the electric field pulse qualitatively follows the current 
pulse. Last, the field from either diffracted Cherenkov or transition is 
proportional to the beam current. Therefore, the radiated power is pro¬ 
portional to the current squared. 
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There are. however, some Important differences between 
Cherenkov and transition radiation. Transition radiation is concerned 
with transiting boundaries. Boundaries associated with electron 
accelerators are the beam port and beam stops. Although data was taken 
with beam stops in position, they will be analyzed in a follow-on 
experiment. The next difference is that Cherenkov radiation occurs in air 
when the beam energy is above the Cherenkov threshold. Also, the 
Cherenkov power radiated (without diffraction) is proportional to the path 
length. Figure 4 [Ref. 4] shows two regions: A and B. Cherenkov radiation 
(without diffraction) is concentrated at one angle (region B). With diffrac¬ 
tion. radiation spreads out about the Cherenkov angle Gc, whereas tran¬ 
sition radiation occurs at all angles (regions A and B) with increased 
intensity in the forward direction. 

B. PREVIOUS EXPERIMENTS 

There have been a number of experiments conducted which measure 
optical, x-ray, or microwave radiation from electron beams. However, rel¬ 
atively few experiments have measured radio frequency radiation from 
propagating electron besuns. Buskirk and Neighbours have previously 
performed a number of such experiments [Ref. 13) dealing with radio 
frequency radiation. The present experiment differs from their earlier 
efforts in the following ways: (1) the number of angles at which data were 
sampled was Increased, (2) many measurements were recorded with 
varied beam parameters, and (3) most importantly, a digitizing camera 
was incorporated into the collection system. This provided a quicker and 


9 









easier method for data collection and storage. This digitizing camera also 
has another important advantage over the previous analog system. The 
resulting data files can be anal 5 rzed by computer programs and softwsure 
without painful conversions. 


C. PURPOSE 

The general overall objective is similar to that of earlier experiments: 
to understand the radio frequency radiation produced an intense elec¬ 
tron beam. More specifically, the purpose of this experiment can be 
defined in the following statements: 

1. Observe differences in propagating to non-propagating beams by 
varying the beam focus. 

2. Attempt to isolate diffracted Cherenkov radiation effect from elec¬ 
tromagnetic pulse (EMP) by varying the beam energy above and 
below the Cherenkov threshold (about 20 MeV). 

3. Obscive possible instabilities in beam propagation due to high cur¬ 
rents. Conversely, investigate the effect low current has on pinching 
the beam. 

4. Analyze the radiation signal and frequency spectrum for empirlral 
obsen'ations. 
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n. EXPERIMENTAL APPROACH 


A. EQUIPMENT ARRANGEMENT 

The specific equipment used in this experiment is listed in Appendix 

A. In order to collect data simultaneously from two antennae and also 
obtain the original source signal, three independent collection systems 
were required. The equipment arrangement for one of the antennae col¬ 
lection systems can be seen in Figure 5. The collection system for the 
source signal is essentially the same, except that one of the accelerator’s 
B dot loops at the beam port replaces the antenna. Actual pictures of the 
collection systems (minus the antennae) are seen in Figure 6. In this 
setup, the radiation from the electron beam is received by the TEM horn 

* antenna and transmitted to the oscilloscope by means of the double- 

shielded transmission cable. The signal on the oscilloscope is digitized by 
a Digitizing Camera System (DCS) and stored on a computer. The B dot 
loop signal was recorded in a similar fashion and was also used to trigger 
the other two oscilloscopes. 

B. ACCELERATOR AND BEAM CHARACTERISTICS 

The accelerator used to produce the electron beam was PHERMEX, 
located at Los Alamos National Laboratory (LAND, New Mexico. A 
schematic of PHERMEX (Figure 7) shows the three cavities which are 
excited by means of radio frequency amplifiers operating at 50 MHz. 

¥ 

PHERMEX can produce a 30 MeV, 900 A beam focused to a 6 mm 
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Figure 7. Schematic of Phermez Accelerator 


diameter. Tlie energy, current, and diameter can be varied as required 
(see Beam Adjustment section). The pulse train created by PHERMEX is 
a maximum of 10 pulses, two to three nanoseconds (ns) long, and sepa¬ 
rated by 20 ns. Pulse trains can be formed up to one per minute [Ref. 
14]. 


C. ANTENNA CHARACTERISTICS 

The transverse electromagnetic (TEM) horn antenna was developed 
at the National Bureau of StandardslRef. 15). The horn is configured as a 
two-conductor, end-fire, traveling-wave structure. Continuously tapered 
reslsUve loading of the two conductive elements provides a relatively 
short, broadband, nondlsperslve antenna which has high directivity. 
Sections of polystyrene foam support the horn’s resistlvcly tapered con¬ 
ductors. The antenna is enclosed In a bakellte box which measures 
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about 30 cm wide by 15 cm high by 51 cm long. A schematic of the 
antenna is shown in Figure 8. The antemia can measure electric fields 
between 1 and 1,600 ^ from 1 MHz to 15 GHz, with a flatness of ±2 dB 
(Ref. 16]. 



Figure 8. Schematic of TEM Horn Antemia 


D. DESCRIPTION OF DIGITIZING CAMERA SYSTEM (DCS) 

The Tektronix DCS (Ref. 17] Is a charge-coupled device (CCD) cam¬ 
era which connects directly to the face of the oscilloscope. 490 vertical by 
384 horizontal picture elements are used to form an image of the trace on 
the oscilloscope. Then the camera transfers the analog video signal to the 
frame store board. The frame store board is a General Purpose Informa¬ 
tion Bus (GPIB) located in an IBM-compatible computer. The conversion 
from analog video signal to digital format is performed by the frame stor¬ 
age board. Here, the digitized raw image which is held in the computer’s 
temporary memory can then be stored on diskettes or a hard disk. The 
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DCS software, along with the frame store board, allows for format con¬ 
version, wave form manipulation, and file transfer. Williams provides a 
more-detailed summary of the DCS S3^tem and software [Ref. 18]. 

E. ANTENNA PLACEMENT 

'Fhe two antennae were each mounted on top of eight-foot ladders. At 
this height, the antennae were approximately two and a half feet higher 
than the beam port. Figure 9 shows one of the actual antenna place¬ 
ments. However, in this picture, both antennae are mounted on the same 
ladder. The positioning of each antenna was relative to the beam path 
from the beam port. At a constant radius of 90 feet, six positions (5-, 10-, 
15-, 30-, 45-, and 60-degree angles) relative to the beam path were mea¬ 
sured and marked. A schematic of the antenna placements is seen in 
Figure 10. Permanent structures on one side of the beam path prevented 
measurements in that direction. 

F. DATA ACQUISITION 

1. Beam Adjustment 

Radiation from the electron beam was recorded in sets of three 
to fiv'e signals per antenna angular position. This procedure was carried 
out to form specific groups comprised of four variations of the electron 
beam (see Table 1). The first group of data sets was taken at full current 
(510 A) and full energy (29.5 MeV). The beam was focused to 6 mm for all 
shots in this group. The second group consisted of data sets received 
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Figure 10. Antenna Placement 


TABLE 1 

DATA COLLECTION 8UMBIARY 



Angles (Degrees) | 

Group 

Beam Paimmetera 

Antenna 

5 

10 

15 

KHl' 

45 

lips!3li 

1 

Full Current & 

A 

X 

X 

X 

X 




Full Energy & 
Tight 

B1 


X 


X 

X 

X 

2 

Half Current & 

A 

X 

X 

X 

1 




Full Energy & 
Tight 

B1 




X 

X 

X 

3 

Full Current & 

A 

X 

X 

X 





Full Energy & 
Diffuse 

B2 





X 

X 

4 

Half Current & 

A 

X 


X 

X 

X 

X 


Half Energy & 
Tight 

B2 

X 

X 

X 

X 

X 
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with the bearn at half current (300 A). All other parameters of the beam 
remained fixed. The third group of data sets was altered only by the 
diameter of the beam. The accelerator's focusing magnet was adjusted to 
deliver a diffused or under-focused beam. In the last group, a half-energy 
(15 MeV), half-current, focused beam produced the signals. 

Throughout all of the shots in this experiment, the pulse train 
remained at either six or seven pulses with a pulse width of 3.5 ns and a 
pulse separation of 20 ns (see Figure 11). Therefore, the pulse trains vai*- 
ied from 121 to 145 ns. 



SIGNAL (Volts) 

Figure 11. Typical B Dot Signal 
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The current of the beam can be determined from the maximum 
peak value. The beam current I is related to the peak voltage V, by 

I * CiV (7) 

where Ci » 129 A/V, V is in volts, and I is in amps [Ref. 14]. 

Current, pulse shape, and number of pulses varied from day to 
day. but consecutive shots were fairly reproducible. 

2. Wavefonn Capture by DCS 

The radiation signal received by the antenna was transmitted to 
the oscilloscope by means of the shielded coaxial cable. The image dis¬ 
played on the oscilloscope was automatically captured by the DCS. The 
signal captured by the camera was immediately digitized by the frame 
store board. The original analog signal was displayed on a TV monitor 
while the digitized signal was graphed on the computer monitor. Prior to 
each Phermex shot, the current oscilloscope settings and intended 
waveform/file name for each waveform were input into the Compaq com¬ 
puter by use of the DCS software. Each waveform was saved on the hard 
disk of the respective computer and a second copy of each file was pre¬ 
served on diskette. 


20 








m. DATA REDUCTION 


A. DATA PILE CONVERSION AND TRANSFER 

The waveforms were recorded and saved In a binary format. Since 
binary Is the normal format when viewing the computer display and 
binary requires less disk storage, this was the preferred format. As previ¬ 
ously mentioned, various waveform manipulations can be accomplished 
using the DCS software In the binary format. However, some waveform 
manipulations such as the Fast Fourier Transform (FFT) cannot be per¬ 
formed with the DCS software. A VAX system was selected to handle the 
additional waveform manipulations. Therefore, an ASCII file format was 
desirable. Nine hundred and seventy-five (975) binary-formatted files 
were individually converted to an ASCII format using the DCS software’s 
built-in feature. During this process, each file was coded with an eight- 
digit file name to provide the maximum information pertaining to the sig¬ 
nal’s capture parameters. 

Once reformatted and coded, all files were transferred from the IBM- 
compatible computer to the VAX system. A KERMIT software program 
provided the means for transferring the files. The VAX system used was a 
Digital Microvax computer 

B. INTERACTIVE DATA LANGUAGE (IDL) PROGRAM APPLICATION 

A program, written with Interactive Data Language (IDL) [Ref. 19] on 
the VAX system, was developed to analyze the collected waveforms. 
Appendix B contains a copy of the main (modified) analysis program. 
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This program corrects for the required In-Une attenuation, the antenna 
response function, and the DC offset. Once the corrections are made, the 
radiation signal and the FFT of the signal are plotted. Directly to the right 
of these plots, tlie corresponding B dot signal and its FFT are depicted. 
Figure 12 is an example of the program’s output. This program was run 
in a batch mode to produce graphs for all the radiation signals. To obtain 
additional information, the program was modified to be user-interactive. 
This interactive program allows the user to truncate the radiation or 
B dot signal to the exact start and finish of the signal. In addition, the 
program will compile the peak frequencies of the FFT and their corre¬ 
sponding relative intensities. 

C. ANTENNA RESPONSE 

Two of the three TEM horns, A and Bl, were calibrated by the 
National Bureau of Standards (NBS). Two sets of responses for each TEM 
horn were tabulated. Each print-out of the data set consisted of 101 
points. The first set contained responses from 50 MHz to 5 GHz at 50 
MHz intervals and a graph of the data (Figure 13). The second set con¬ 
sisted of responses from 1 MHz to 101 MHz at 1 MHz intervals and a 
graph of the responses from 1 MHz to 512 MHz (Figure 14). The overall 
uncertainty of the calibration is ±1.0 dB. 

A print-out with the response for eveiy MHz from 1 MHz to 600 MHz 
would have been ideal. However, an attempt was made using the second 
set’s graph (Figure 14) to extract the response from 5 MHz to 512 MHz at 
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RELATIVE INTENSITY SIGNAL (Volts) 



Figure 12. Example of Program Output 
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Frequency (MHz) 
Figure 14. Anteniw Response 
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5.3 MHz intervals. This was referred to as the “fine correction." The val¬ 
ues provided by the calibration print-out corresponding to Figure 9 were 
referred to as the “coarse correction." 

The c!orrection values are referred to by the NBS National Engineer¬ 
ing Laboratory (Ref. 201 as the TEM horn transfer function. The transfer 
function is defined as 


TF = 20 LOG g 


( 8 ) 


where 

V = output voltage of TEM horn, rms volts, and 
E = electric field strength, volts/meter 
The transfer function correction using equation 8 was included in 
the analysis program to calculate the irTTs. Points lying between the 5 
MHz Interval in the fine correction case and between 50 MHz interval in 
the coarse correction case were derived by linear Interpolation. 

A comparison of the fine correction to the coarse correction is seen 
in Figure 15. The uncorrected FFT is also plotted to show the effect of the 
transfer function. The uncorrected FFT displays the characteristic 50 
MHz peak and its harmonics. The coarse correction has a similar trace 
but it was Increased in Intensity by a factor of about 100. However, the 
fine correction causes a flattening of the characteristic peaks. This effect 
is attilbuted to rapid fluctuations in the graph of responses and to the 
method of evaluating the closest value. The fine correction had an added 
disadvantage: it lengthened the program so much that it could not be 
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executed in IDL. The coarse correction was therefore used in evaluating 
all the FFTs. Figure 16 depicts the transfer function for both antenna A 
and antenna Bl. 

Also from the MBS test procedure, the voltage received is related to 
the electric field E at the antenna, by 

E « CfV (9) 

where Cf = 96.6, V is in volts, and E is in volts/m [Ref. 201. 
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Figure 16. Antenna A & B1 Transfer Functions 
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D. SIGNAL REPRODUClBILnT 

During the experiment, the reproducibilit]/^ of signals was investi¬ 
gated. A set of three signals obtained at the same angle, under the same 
beam parameters, and with the same antenna (A) are shown in Figure 
17. One immediately notices that the three signals are very similar. How¬ 
ever, minor but observable differences in the peaks are apparent. Tlie dif¬ 
ferences can be attributed to slight changes in the electron beam coupled 
with an uncalibrated DCS system and slight DC offsets. In general, the 
recorded signals for each data set are reproducible. 

Figure 18 is a set of three signals taken under the same conditions 
as Figure 17, except that antenna B1 was used. Here, similar effects are 
seen. Again, the overall conclusion is that the recorded signals are repro¬ 
ducible in each data set. 

E. ANTENNA-TO-ANTENNA COMPARISON 

As seen in Table 1, antenna A was used primarily at the smaller 
angles of 5, 10, and 15 degrees. B1 or B2 was normally placed at the 
larger angles of 30, 45, and 60 degrees. Tlierefore, a comparison of 
antenna A to antenna B was required to relate the smaller- and larger- 
angled values. Figure 19 depicts the signals from antennae A and Bl. 
The beam shot and angular position were exactly the same for both 
antennae. The signals display similar oscillations and characteristics. 
However, the maximum peak-to-peak voltage of antenna Bl is approxi¬ 
mately one-fifth of antenna A. Antenna B2 was also compared to Antenna 
A in a similar manner. Table 2 shows the antenna conversions for Bl 
and B2 to antenna A. 
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Figure 19. Comparison of Antenna A and B1 Signals 


TABLE 2 

ANTENNA CONVERSION 


Antenna 

Factor to Convert Antenna B 
Signal to Antenna A 

CPC 101 (Antenna Bl) 

4.9 

CPC 102 (Antenna A) 

1.0 

CPC 103 (Antenna B2) 

0.33 


It was found that antenna B1 compared much more favorably than 
antenna B2 to antenna A. The signals from antenna B2 were very sharp 
and appeared almof^t erratic. Antennae A and B1 signals were much 
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smoother. However, as mentioned in the previous section, the signals 
^vere reproducible for each data set recorded v/lth the same antenna. 


31 










IV. RESULTS ANt> DISCUSSION 

' I ■ . 

A. RESULTS 

The experimental results are organized into five areas: signals and 
related spectrum, signal strength as a function of angle, relative intensity 
of spectrum harmonics, shifting of fundamental frequency, and modula¬ 
tion of spectrum. A discussion of tlie following results appears in the 
next section, 

1. Signals and Related Spectrum 

A complete set of all radiation and B dot signals along with their 
respective spectra was created using the IDL program. However, here one 
signal per angle (5. 10, 15. SO, 45. arid 60 degrees) and beam parameter 
are provided for completeness. Figure 20-22, Figures 23-25, Figures 26- 
28. and Figures 29-^1 are signal and B dot results for full current, half 
current, diffuse beam (full current), and half energy beam (half current), 
respectively. 

2. Signal Strength as a Function of Angle 

For each radiation signal from each data set, the maximum 
peak-to-peak voltage Avas measured. This value divided by two was also 
converted into the actual electric field strength using equation 9. The 
current for each beam shot was computed using equation 7 and tabu¬ 
lated in Table 3. Difficulties with the B dots on the first day caused 
inaccurate current readings. It Is assumed that the current was 
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Figure 20. Signal and Spectrum for Full Current, Pull Energp, 
Angles of 5 and 10 Degrees 
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Figure 22. Signal and Speotnim for Pull Current, Full Energy, 
Anglea of 45 and 60 Degrees 
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Figure 23. Sigaal and Spectrum for Half Current, Full Energy, 
Angles of 5 and 10 Degrees 
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Figure 27. Sigiud find Spactmm for Dlffote Beam, Full Current, 
Full Bnerg 7 , Angles of 15 and 30 Degrees 
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Figure 28. Signal and Spectmm for DilAiae Beam, Full Current, 
Pull Energy, Anglea of 45 and 60 Degrees 
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Figure 31. Si^ial and Spectnim for Half Energy, Half Ctxrrent, 
Angles of 45 and 60 Degrees 
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TABLES 

ACTUAL BEAM CURRENTS ON AMPS) 


Beam 

Parameter 

ABflea hi Degreea | 

10 

10 

IB 

HEOI 

45 

60; ) 

Full Current 

495 ± 13 

286±12* 

275 ±. i4‘ 

265 ± 5‘ 

270 i. 2* 


Half Current 

302 ± 15 

404 ±2 

297 ±2 

265 ± 5 

610 ±7 

292 ± 13 

DifTuse 

463 ± 22 

449 ± 5 

456 ±5 

463 ± 22 

449 ±5 

466 ±5 

Half Energy 

312 ±2 

288 ± 13 

268 ± 13 

303 ± 20 

319 ± 1 

303 ± 7 


‘Problems with B dot registered lii accurate currents. Actual currents were approxi¬ 
mately 500 amps. 


approximately 500 A for that day (for full-current tight beam). Since the 
current, angle, and distance are ail known values, the expected electric 
field can be determined using equation 6. A comparison of expected to 
actual electric fields is in Table 4. The average radiation peak-to-peak 
value was also plotted as a function of angle. Full versus half current 
(Figure 32), tight versus diffuse beam (Figure 331, and full versus half 
current (Figure 34) were compared. In addition, the ratio of the upper 
curve to the lower curve for the plotted points was compared. The ratio of 
curves for full versus half current, tight versus diffuse beam, and full 
versus half energy are exhibited in Figures 35, 36, and 37, respectively. 

3. Relative Latensity of Spectnun Harmonici 

The fundamental frequency and harmonics are easily observed 
in the spectrum. By using the interactive version of the IDL program, the 
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TABLE 4 

COMPARISON OF ACTUAL TO SXPRCTSD ELECTRIC FIELD 



Degvses 

S 

10 

18 

30 

48 


Pull 

Cunvnt 

4 

Actual 

mj 

U.4 ± .7 

12 ±1 

4.91.6 

2.01.3 

2.121.06 

2.3 .-..1 


Ejected 

^ mj 

9.85 ± .09 

6.41.1 

4.21.1 

2.081.02 

1.351.01 

1.001.01 

Half 

Current 


Actual 

< mj 

9.85 ± .09 

3.31.1 

2.61.1 

B 

1.51.1 

1.32 :b .05 


Eiepected 

H 

5.041.02 

2.461.02 

1.271.03 

.82 1.02 

.56 1 .o:t 


relative intensity for each harmonic J was then plotted versus angle:. Thi s 
was repeated for the first four higher harmonics (Figures 38-41) for the^ 
four different beam parameters. However, the full current was carried ouv. 
to the ninth harmonic (see Figure 42). 

4. Harmonic Frequency Shift 

Again, by using the interactive version of the IDL program, the; 
frequency of eacli harmonic was evaluated. The average harmonic for 
each shot was calculated by plotting the frequency against the harmonic 
number J. A least square fit of the points provided the slope or the 
overall average frequency. The slopes for shots at the same angle and 
beam parameters were averaged. This average frequency (approxlmsitely 
50 MHz) was then plotted versus angle for both full (Figure 43) and half 
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Figfire 32. Signal Strength of Full Versue Half Current (Pull Energy) 
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Figure 33. Signal Strength of Tight Versus Difhise Beam 
(Full Current* FuU Energy) 
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Figure 34. Signal Strength of Pull Versua Half Energy (Half Current) 
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Figure 35. Ratio of Full to Half Current Signals 
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Figure 36. Ratio of Tight to Diffuse Beam Signals 
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Figure 37. Ratio of Full to Half Energy Beam Signals 
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Figure 42. 5«0 Karmonict of Full Current Spectrum (Full Energy) 
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Figure 43. Average Frequency of Spectrum for Full Current 

Figure 44) current. The same calculations were performed on the B clot 
spectrum, pnd the respective plots are in Figure 45 and 46. To determine 
whether shifting of the harmonic frequency had occurred, the B dot full 
current average frequency (Figure 45) was subtracted from the full cur¬ 
rent average frequency (see Figure 47). Similarly, the difference was eval¬ 
uated and plotted for the half current values (Figure 48). 

5. Modulation of Spectrum 

The FFT results of the B dot are similar to the results of the pat¬ 
tern seen from a diffraction grating in optics. From optics [Ref. 21], the 
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Figure 46. Average Frequency of Spectrum for Half Current B Dot 
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Figure 47. Signal and B Dot Difference for Full Current 









where c is the speed of light and d is the slit width or pulse width. A clear 
example of modulation is observed in Figure 49. Here, a simulated B dot 
pulse is formed by seven equal-volt 2.6 ns pulses separated 20 ns. In 
this case, the pulse width (2.6 ns) acts as the slit width and is related to 
the modulation minima. The corresponding minimum is therefore the 
reciprocal of 2.6 ns, which is 385 MB^. The 20 ns pulse separation forms 
the 50 MHz peaks and higher harmonics. 
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Figure 49. Modulation Example from a Perfectly Pulaing Accelerator 


Modulation minima were also observed In the actual spectrum of the 
signal and B dot. An actual B dot and its spectrum are seen in the top 
right section of Figure 24. Although the beam is not perfect, it exhibits 
similar modulation minima at a slightly lower frequency (360 MHz). Table 
5 lists first, second, and third modulation minima for various angles and 
beam parameters. These values have up to a five percent margin of error. 
In addition, only those signals listed in Table 5 had clear enough modu< 
lation minima co discern a value. 
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TABLES 

MODULATION MINIMA 


Beani 

Parameter 


1st 

2nd 

ard 

Full 

Current 

BDct 

275 

525 

750 

5® 

270 

450 

610 

B Dot 

275 

580 

905 

10® 

280 

— 

— 

B Dot 

285 

595 

890 

60® 

290 

600 

895 

Half 

Current 

BDot 

360 

705 

895 

15® 

320 

675 

— 

B Dot 

275 

500 

720 

60® 

295 

600 

895 

Half 

Energy 

BDot 

220 

470 

740 

5® 

220 

475 

720 

BDot 

220 

470 

740 

10® 

240 

450 

650 

BDot 

230 

505 

795 

15® 

235 

495 

755 

B Dot 

— 

445 

860 

30® 

— 

440 

815 

B Dot 

— 

445 

795 

60® 

— 

450 

740 


B. DISCUSSION 

The five areas outlined In the results section will now be discussed. 
Some areas overlap and are inserted wherever clarification is necessary. 
1. Signals and Related Speetnun 

There are a number of general observations which hold for all 
signals independent of antenna angle or beam parameter. First, the 
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overall pulse train length of the radiation signal and the B dot signal are 
essentially equal. Second, since the radiation depends on the rate of 
change of the current, there is a one-to-one correspondence of B dot 

, peaks to major radiation signal oscillations. One can also observe the 

successive differences in B dot pulse peaks as a similar difference in the 
successive peak oscillations. Last, the 50 MKz frequency and higher 
harmonics are clearly evident in all radiation spectra. However, the dif¬ 
fuse beam radiation spectrum harmonics were less consistent and dis¬ 
tinguishable. Also, the modulation minima of the diffuse beam were not 
clearly evident from the spectrum. 

2. Signal Strength a* a Function of Angle 

By viewing Figures 32. 33. and 34, one discovers a general 
decrease in signal strength with increase in angle. This effect is mainly 

w 

attributed to transition radiation. A comparison of the general transition 

* radiation pattern (Figure 3) and the radiation (Figures 32-34) observed 
from this experiment shows significant similarities. Actual electric field 
strength from Table 4 is comparable to the expected strength at angles of 
15 degrees or greater. However, at smaller angles, the actual electric field 
is considerably larger. This increase in field strength may be attributed to 
diffracted Cherenkov radiation and diffracted transition radiation. A 
closer look at the 5- and 10-degree values of both curves on Figure 34 
shows a large difference at these angles. The full energy (29.5 MeV) is 
above the energy to cause Cherenkov radiation, and the half energy (15 
MeV) is below the required energy. The ratio from Figure 37 at small 

• angles also suggests diffracted Cherenkov radiation may be present. 
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The ratio of the full current peak-to-peak signal versus the half 
current peak>to-pcak signal (Figure 35) exhibits (with the exception of the 
10-degree data point) a fairly constant ratio of approximately 1.4. On re¬ 
examining equation 6. one finds that the ratio of full versus half current 
signals at the same angle is merely the ratio of the currents. All other 
values in equation 6 remain constant and, therefore, from Table 5 the 
average ratio of full to half currents is approximately 1.6. 

In Figure 33. a propagating (tight) beam and a non-propagating 
(diffuse) beam are compared. Again, one observes a similar transition 
radiation pattern. lx>oklng again at Figure 36, one notices a small ratio 
(of approximately 2) for the smaller angles and a larger ratio (of approxi¬ 
mately 8) for the larger angles. It appears that a non-propagating beam 
still has significant radiation in the forward direction but radiates signifi¬ 
cantly less than a propagating beam in the other directions. 

3. Relative Intensity of Spectrum Harmonics 

As previously mentioned, the 50 MHz harmonics are readily 
observable in the radiation spectrum. Figure 38-41 indicate the relative 
intensity of 50 MHz (J » 1) and the next three higher harmonics. In gen¬ 
eral, the harmonic intensity follows the same pattern as the signal 
strength. Therefore, increases in angle cause a decrease in relative 
intensity. Also, each successive harmonic curve decreases in intensity. 
The relative intensity is still significant in higher order (sometimes up to 
J B 20) harmonics. This effect is seen in each signal FFT and is graphi¬ 
cally depicted in Figure 42 with J ■ 5 through 9. 
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4. Hftfmonio Frequency Shift 

Figures 45 and 46 indicate that a relatively constant 50 MHz 
beam was being produced. Viewing Figures 43 and 44, a frequency shift 
with increasing angle is Indicated. To negate any effect changes in the 
beam source had on this shifting, the average B dot frequency was sub* 
tracted from the average signal spectrum frequency {Figures 46 and 47). 
Again, a slight shift to higher frequency with increasing angle is evident. 
This predicted pattern was seen earlier in Figure 2. These experimental 
findings seem to support the theory behind Figure 2. However, measure¬ 
ments In both Figures 43 and 44 Involve two different collection systems. 
One collection system captured the signals at 5, 10, and 15 degrees, 
whereas the other system collected the signals at the remaining angles. 

5. Modulation of Spectrum 

The summary of modulated minima (Table 5) excludes minima 
from the diffuse beam. As mentioned previously, modulation minima 
were difficult to locate for diffuse beams and were consequently omitted. 
At larger angles, the modulation minima were easier to discern. This fact 
is in agreement with the BuskirK, Neighbours, and Maruyama model 
[Ref. 9]. From the remaining beam parameters, there is significant 
agreement between the B dot and signal spectrum minima. In general, 
the first and sometimes the second minima are consistent. The first min¬ 
ima appears on the average at 235 MHz for a five-pulse beam, 280 MHz 
for a six-pulse beam, and 340 MHz for a seven-pulse beam. Using equa¬ 
tion 10 and multiplying by the respective number of beam pulses 
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consistently gives values between 3.4 and 3.6 ns. Hiese values are very 
close to the actual value of the beam pulse width (3.5 ns). 
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V. CONCLUSIONS 


Radio frequency radiation, produced by the electron beam of the 
PHERMEX accelerator, consisted mainly of transition radiation (EMP) 
and minor contributions from diffracted Cherenkov radiation (when at 
full energy). Regardless of the beam energy, current, or frequency, the 
following general observations can be made concerning the analyzed 
radiation and its associated spectrum: 

1. There is a one-to-one relation between the radiation oscillating 
peaks and the electron beam pulses. 

2. The length of the radiation pulse train remained constant and con¬ 
sistent with the beam pulse train. 

3. PHERMEX operating frequency of 50 MHz was consistently 
observed with up to the 20th harmonic in the signal spectrum. 

4. In addition, there is some evidence that shifting may occur to the 
fundamental frequency and its harmonics in the signal spectrum. 
Increasing angle seems to cause an increase in the average har¬ 
monic. 

5. Electric fteld strength of the radiation is proportional to the beam 
current. 

6. Modulation minima of radiation signal spectrum were similar to B 
dot spectrum. The first modulation minima is directly related to the 
beam pulse width. 

In addition, differences between propagating and non-propagating 
beams were noted as follows: 

1. Propagating beams retain significant beam characteristics through 
their spectrum, whereas non-propagating beams do not. 

2. Non-propagating beams produce less radiation in the forward direc¬ 
tion than a propagating beam. 
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3. A signlilcant decrease in radiation exists between propagating and 
non-propagating beams in other than the far forward direction. 

There has been a major improvement in the collection system. Sig¬ 
nals can now be quickly digitized and stored in a computer for later anal¬ 
ysis. Previous experiments resulted in only analog measurements. 

Further work, both theoretical and experimental, is already planned 
to continue investigating this and additional data collected during the 
PHERMEX experiment. 
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APPENDIX A 

EQUIPBISNT LIST 


1. Tektronix 7103 Oscilloscope (3 each) 

2. Transverse Electromagnetic Horn Antenna (3 each) 

3. Tektronix Digitizing Camera S3rstem (DCS) (3 each) 

4. Personal Computers 

Compaq Portable II Model 4 (2 each) 

IBM XT (1 each) 

5. Video Monitors (3 each) 

6. Double-Shielded Cables (3 each) 

7. Tektronix Video Processor HCOZ 

8. IBM Graphics Printer 
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APPENDIX B 

IDL ANALYSIS PROGRAM 


; T!>ia is an updatad varaien at rf_axp'c by Rich Lally on 9 APR 90. 

,'PROSRAM TO AHALXZE RT DATA TAXBH MITH DCS 
FOR Z«1,S DO BBGZN 
PRINT,' ' 

CNOrOR 

PRINT," PROORAM PLOTTER.PRO" /Profraai Tltla 

PRINT," " 

PRINT,"*••• Pre^rac to analyaa DCS Data. ***** 

FOR I-1,S DO BEGIN ;Skips 5 linaa 

PRINT,' ' 

ENDFOR 

READ,'DO YOU NISH TO FIND MAXIMUMS 7 1 FOR TBS AND 2 FOR NO',FMAX 


PRINT,"Fila MUST NOT hava an .ASC axtansaon." 

immmmmmmmmmmm INPUT OF FILENAME mmmmmmmmmmmmmm m mmmmmmmmmmmmmm m mmmmmmmm 

NAME-' ' 

R£AS,'Entar tha naaa o£ tha data fila',NAME ,'Raaata fila naaa to null 

; Sat 

FOR JJ-l,2 DO BEGIN 


;Maaaursa langtb of fila nana 


/INPUTS CORBSPONDZNG BDT NAME 
;8tripa off disaetory t .ase 

DATLETTERoSTRMIO(TITLE,0,1) 

CASE DATLETTER OF * 

'M': DAT-'HON' 

'T': DAT-'TUB' 

'N': DAT-'NED' 

'H'I DAT-'THU' 

'F'! OAT-'FRI' 

ELSE: DAT-'UNRNONN' 

ENDCA8E 


NAMBLBN-8TRLEN(NAME) 

IF JJ EQ 2 THEN STRtUT,NAME,'BDT',1 
TITLE-STRMID(NAME,0,(HAMBLEN)) 


DEORXBS-STRMIO(TITLE,1,2) 

ANTIMNA-STRMZD(TITLE,3, X) 

SHOTireMBBR-ETRMIO(TITLE. 4,3) 

CORREHTL-STRMID(TITLE,7,1) 

CAM CURNDITL OF 
'F'! CORKENT-'FULL' 

'B': CORMBIT-'IIALF' 

ELSE: COMUtNT-'UNlOIONN' 

SmCASE 


IF OEeiVEBE EQ 'BO' TEEN TOPLZNB-TITLE a * ODRRMT * ' I,' S 
* 'B DOT FM '*DAT SHOTNUNEER ELBE f 

TOPLZNE-TZTLB * '•' * CUMOENT *■ ' Z,' * DEGREES * ' DEG,AMT ' t 
* ANTEMNA * ',' *DAY * ' 9' * SKOTNUMBER 

FZLB - KAIB * ' .ASC' 

;-m. . saiiaiiiis READ DATA IN FILE ■ ■■■as a s a iiaa — ■ 
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;This is an updated version of rf_«xpt by Rich Laliy on 8 APR 90. 

; PROGRAM TO AMAZ.TZE RF DATA TAKEN WITH DCS 
rOR 1-1,5 DO BEGIN 
PRINT,' ' 

ENDFOR 

PRINT," PROGRAM PLOTTER.PRO" ;Program Title 

PRINT," " 

PRINT,Program to analyze DCS Data. ****" 

FOR 1-1,5 DO BEGIN /Skips 5 lines 

PRINT,' ' 

ENDFOR 

READ,'DO TOO WISH TO FIND MAXIMUMS ? 1 FOR YES AND 2 FOR NO',FMA3C 


PRINT,"File MUST NOT have an .ASC extension.” 

I mmmmmmmmmmmm INPUT OF FILENAME mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmw^ 

NAME-'' 

READ,'Enter the name of the data file',NAME /Resets file name to null 

/Set 

FOR JJ-1,2 DO BEGIN 


nAMELEN-STRLEN(NAME) /Measures length of file name 


IF JJ EQ 2 THEN STRPUT,NAME,'BDT',1 /INPUTS CORESPONDING BDT NAME 

TITLE-STRMID(NAME,0,(NAMELEN)) /Strips off directory & .asc 

D,.iLSTTER-STBMID (TITLE, 0,1) 

CASE DAYLETTER OF 
'M': DAY-'MON' 

'T'! DAY-'TUE' 

'M': DAY-'WED' 

'H': DAY-'THU' 

'F': DAY-»'rRt' 

ELSE 1 DAY-'UNKNOWN' 

ENDCASE 


DEGREES-STRMIO(TITLE,1,2) 

AMTENNA-STRNID(TITLE.3,1) 

SHOTNUMBER-8TRMID(TITLE,4,3) 

CURRENTL-STRMID(TITLE,7,1) 

CASE CURRENTL OF 
'F'; CURRENT'-'FULL' 

'H': CURRENT-'HALF' 

ELSE: CURRENT-'UNKNOWN' 

ENDCASE 


IF DECREES EQ 'BD' THEN TOPLlNE-TlTLE + '-' + CURRENT + ' I,' S 
* 'B DOT FOR '•t’DAY 4- SHOTNUMBER ELSE $ 

TOPLINE-TZTLS 4- CURRENT * ' I,' 4- DEGREES * ' DEG, ANT ' $ 

* ANTENNA 4- ', ' 4-DAY 4- ' »' 4- SHOTNUMBER 

FILE - NAME 4- ' .ASC' 
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0PENlR,2,FIliE ;Opmnt input fil« for rtading 

hoadorvstrsrrieOfSO) /Oinanaidnn rtrlng «rr«yi 

rdhdroatrarr( 80 , 1 ) 

DATA « F]:.TAIW(512) 

AA " ' ' 

;R*acl in valuas from haadar 


I-O 

tmiLE NOT EOF (2) DO BEGIN 

IF ((I GT 19) AND (I LT 532)) THEN BEGIN 
READF,2, '$(Fe.6) FF 
DATA (1--20) « FF 
SNDXF ELSE BEGIN 
REAOF,2, AA 
KEADER(I) - AA 
ENDELSE 
I-I + 1 
ENDWHILE 


':') +2 




(HEADER(6) , 

POSl, 

3) , 

2) 

':') +2 




(HEADER(7), 

POS2, 

3), 

2) 


POSl •• STRFOS (KEAOEB(6), 

STARLBL - STRTRIM (STRMID 
POS2 • STBFOS (HEADER(7), 

STOPLBL - STRTRXM (STRMXD 
START - FIX(STARLBL) 6 STOP • FIX(STOPLBL) 
DELTA-STOP-START + 1 
POS3 « STRPOS (HEADER(12), ';') +2 
TXKLBL « STRTRXM (STIRMID (HEADER(12) ,POS3,20) , 2) 
TIME - FLOAT(TIMLBL) 

IF (DEGREES EQ 'BO'} AND (TIME GT S.S-9) THEN TIME 


STRING TO NUMBER CONVERSION 


4.34783B-10 


T-incl 9 an(<lalta) *tiaa /TIME DATA 

; ■—w—^.■■■■■ m» «attsnuatiow correction — u— 

OB-0 

IF (DATLBTTER EQ 'M'l AND (ANTENNA EQ 'A') THEN DB-100.0 

IF (DATLETTER NE 'M') AND (ANTENNA EQ 'A'} THEN DB-63.1 

IF (OATLSTTER EQ 'M') AND (ANTENNA EQ 'B') THEN OB-3.2 

IF (OATLBTTEK EQ 'T') AND (ANTENNA EQ 'B'> THEM DB-IC.O 

IF (DATLBTTER EQ 'W') AMD (ANTENNA EQ 'B') THEN DB-10.0 

IF (OAYLETTER EQ 'F') AND (ANTENNA EQ 'B') THEN DB-316.2 

IF (DATLETTER EQ 'H') AND (ANTENNA EQ 'B') THEM BEGIN 

IF (SKOTNUMBER LT 019) THEN DB-10.0 

IF (£HOTNUNBCR GE 019) AND (SHOniUHBER LT 027) THEN OB-100.0 
IF (SKOTNUMBER GE 027) THEN DB-316.2 
ENDIF ELSE BEGIN 

IF (DEGRESS EQ 'BO') THEN DB-70.711 

ENDELSE 


cloaa,2 


Ixtitla-'TIMS (naae)' 

IYtitla-'SIGNAL (Volta)' 

Itypa- 16 
IMTITLE-TOFLINE 

X-T*1.0B9 /SCALE TO NANOSECONDS 

UY-DATA(START:STOP) 

Y-UY * OB /ATTENUATION CORRECTION 

/mmrnmrnmmmmmmOJKSOR ROUTINE TO SISB WAVEFORM — 

SET PLOT,l 

IF ?J EQ 2 THd GOTO,LIiT 
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!C-0 


PL0T,X,5r 

GOTO,iirSR 

LIFE: LX«X & LV-Y 

ERASE 

TIME (nice)' 

!ytitl«o'SIGNAL (Volta)' 
!typ«» 16 
IMTITLE-TOPLII^ 
!XMAX-MAX(X) 

!C-0 
!GRID-0 
lXTICKS-0 


PLOT,LX,LY 

LIFER;SEEEEEE-0 

CURSOR,SI,S2,XI,Y1 

CURSOR,S3,S4,X2,Y2 

CORRECT-TIME*!.0E09 

X2-X2/CORRECT 

Xl-Xl/CORRECT 

IF XI LE 0 THEM XI - 0.0 

XXX-X(XI:X2) & YYY-Y (XI:X2) 

DCOFFSET»Y(Xl) 

Y-YYY-DCOFFSET 
ZEROER-X(XI) 

X-XXX~ZEROER 
OELTA-N_ELEMENTS(X)-1 

; ————————— FOURIER TRANSFORM ■■■■■■■■■■■wwiiim— 

I-IMDGEK(DELTA) 

MU-1/(TIME*DELTA) /FUNDAMENTAL FREQUENCY 

P2-OELTA/2 

FREQ-FLTARR(P2) 

CDB«FLTARR(P2) 


FOR J-1,P2-1 DO BEGIN 

FREQ(J)-I(J)*NU/1.0E6 ;IN MHZ 

BAO-HHSRB ( FREQ LE 600 ) ;TO CHANGE MAX FREQ CHANGE NUMBER I.E. 600 
BAD-WKERE ( FREQ LE 1200 ) ;TO CHANGE MAX FREQ CHANGE NUMBER I.E. 600 

'————ANTENNA CORRECTION——————————— 

IF DEGREES EQ 'BD' THEN GOTO,DIVE 
IF ANTENNA EQ 'A'THEM BEGIN 


CASE 1 OF 

(FREQ(J) GE 1) AND (FREQ(J) LT-SO): CDB(J)-((FRBQ(J)-1)/(50- $ 

1))*(30.1-30.0)+30.0 

(FREQ(J) GB 50) AND (FREQ (J) LT 100): CDB (J) - ((ITRBQ (J)-50) / (100- $ 

50))*(30.3-30.1)+30,l 

(rREQ(J) GE 100) AND (FRaQ(J) LT 150): CDB(J)-((FREQ(J)-100)/(150- 
100))*(37.9-30.3)+30.3 

(FREQ(J) GE 150) AND (FREQ(J) LT 200): CDB(J)-((FREQ(J)-150)/(200- 
150))*(41.9-37.9)+37.9 


(rRSQ(J) GE 200) AND (FREQ(7) 
200))*(32.1-41.9)+41.9 
(FRXQ(J) GE 250) AND (FRBQ(J) 
300))*(30.7-32.1)+32.1 
(FREQ(J) GE 300) AND (FREQ(J) 
300))*(36.0-30.7)+30.7 
(FREQ(J) GB 350) AND (FREQ(J) 
350))*(34.8-36.0)+3S.0 
(FRSQ(J) GB 400) AND (FRBQ(J) 


LT 250): CD8(J)-((FREQ (J)-200)/(250- 
LT 300): CDB(J)-((FREQ(U)-250)/(300- 
LT 350): CD8(J)-((FREQ(J)-300)/(350- 
LT 400): CDB(J)-((FREQ(J)-350)/(400- 
LT 450): CDB(J)-((FREQ(J)-400)/(4S0- 


$ 

$ 

$ 

$ 

$ 

$ 

$ 
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400) ) * <37.3*-34.8)+34.8 

(riUCQCJ) CB 450) AND (FNSQ(J) LT SCO) : CD8 ( J) - ((flUCQ'Cr) -450) / (500> $ 

450))*(38.5-37.3)+37.3 

(FNEQ(J) GE 500) AND (FnBQ(J) LY 550): CDB C J) «<(FMBQ( J) -500)/(550- $ 

500))*(42.5-38.5)+38.5 

(FRBQ(J) GE 550) AND (FREQ(J) LT 600): COB (J)•((FlUBQ(J)-550)/(600- $ 

S50))*(39.6-42.5)-t.42.S 


ELSE: CDB(J)-40.0 
ENDCASE 

ENDXF ELSE BEGIN 

IF (DAY EQ 'H') OR (DAY EQ 'F') THEN GOTO, DIMER 
IF (DAY EQ 'W') AND (SMOTNUMBBR 6S 008) THEN GOTO,DXMBR 


CASE 1 OF 

(FREQ(J) GE 1) AND (FBEQ(J) LT 50): CDB(J)-((FREQ(J)-1)/(50- S 
X))*(30.7-33.3)+33.3 

(FREQ(J) GE 50) AND (FREQ(J) LT 100): CDB(J)-((FREQ(J)-50)/(100- $ 

50))’*(33.4-30.7)+30.7 


<FREQ(J) GE 100) AND (FREQ(J} LT 
100)) * (54.0-33.4)+33.4 
(FREQ(J) GE 150) AND (FREQ(J) LT 
1S0))*(34.6-54.0)434.0 
(FREO(J) GE 200) AND (FREQ(J) LT 
200))«(32.8-34.6)434.6 
(FR£Q(J) GE 250) AND (FREQ(J) LT 
300))*(35.0-32.8)432.8 
(FREQ(J) GE 300) AND (rREQ(J) LT 
300))*(37.5-35.0)435.0 
(FREQ(J) GE 350) AND (FREQ(J) LT 
350))*(35.3-37.5)437.5 
(FREQ(J) GE 400) AND (FREQ(J) LT 
400)) <*(37.4-35.3)435.3 
(FRSQ(J) GE 450) AND (FREQ(J) LT 
450))'*(38.4-37.4)437.4 
(FREQ(J) GE 500) AND (FRBQ(J) LT 
500))*(42.3-38.4)438.4 
(FREQ(J) GE 550) AND (FREQ(J) LT 
550))'*(41.7-42.3)442.3 


150): CDB(J)«((FREQ(J)-100)/(150- 
200): CDB(J)-((FREQ(J)-150)/(200- 
250) : CDB (v-*)-( (FREQ (J)-200)/(250- 
300): CDB(J)-((FREQ(J)-250)/(300- 
350): CDB(J)-((FREQ(J)-300)/(350- 
400): CDB(J)-((FREQ(J)-350)/(400- 
450): CDB(J)-((FREQ(J)-400)/(450- 
500): CDB(J)«( (FREQ(J)-450)/(500- 
550): CDB(J)-((FREQ(J)-500)/(550- 
600): CDB(J)-((FREQ(J)-550)/(600- 


$ 

5 

6 
s 
$ 
$ 
$ 
$ 
$ 
$ 


ELSE: CDB(J)>«40.0 
ENDCASE 


ENDELSE 
GOTO,FLIP 
DXMFR: CASE 1 OF 

(FRBQ(J) GE 1) AND (FREQ(J) LT 50): CDB(J)-((FREQ(J)-1)/(50- $ 

1))4(51.0-31.5)431,5 

(FRF.Q(J) GE 50) AND (FREQ(O') LT 90): CDB (J) > ((FREQ (J)-50) / (90- $ 

50))*(44.5-51.0)451.0 

(FRBQ(J) GE 90) AND (FREQ(Ji LT 100): CDB(J)-((FREQ(J)-90)/(100- $ 

90))*(48.0-44.5)444.5 

(FREQ(J) GE 100) AND (FREQ(J) LT 200) : CDB (J)-((FREQ (J)-UO) / (200- $ 

100))*(55.5-48.0)448.0 

(FREQ(J) GE 200) AND (FRBQ(J) LT 300): CDB(J)«((FREQ (J)-200)/(300- $ 

200))4(35.7-55.5)455.5 

(FRBQ(J) GE 300) AND (FREQ(J) LT 400): CDB(J)-((FREQ(J)-300)/(400- i 
300)) *(52.0-55.7)455.7 

(FRI,Q(J) GE 400) AND (FRBQ(J) LT 500) : CDB(O) ■( (FREQ (J)-400) / (500- $ 

400))*(47.0-52.0)452.0 

(FREQ(J) GE 500) AND (FREQ(J) LT 600); CDB(J)-((FREQ(J)-500)/(400- $ 

500))*(41.5-47.0)447.0 

(FREQ(J) GE 600) AND (FRBQ(J> LT 700): CDB(J)«((FREQ(J)-600)/(700- $ 

600))*(48.0-41.,•5)441.5 

(FRBQ(J) GE 700) AND (FRBQ(J) LT 800); COB(J)•((FREQ(JJ-700)/(800- $ 
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700))• (49.0-4e.0)+4e.0 

(FIUEQ(J) GS 800) AND (FRSQ(J) tT 900); CDB(J) •( (rNBQ(J) -800) / (900- $ 

800) )* (43.5-49.0)4-49.0 

(FREQtJ) GS 900) AND (FASQ(.7) LT 1000): COB (J)-((nUBQ (J)-900) / (1000- S 
900))»(43.S-43.5)+43.5 

BL8K: COB(J)-43.5 


SNDCASr 

GOTO,FLIP 
DIVE: COB(J)-0.0 
FLIP: FFF-1.0 
SNDFOR 


FBEQ(0)«0.0 

IF ANTBltlNA EQ 'A'THEN COB (0)-30.0 
IF ANTENNA BQ 'B'THEN COB (0)-33.3 
BADZE-MAX(BAD) 

;WITHOUT ANTENNA RESPONSE 
;ANTENNA RESPONSE FUNCTION 


;Z-FFT(Y,1) 
CT-10'' (CDB/20) 


02-FFT (Y,1) 
Z«U2*CT 


/WITH ANTENNA RESPONSE 
/WITH ANTENNA RESPONSE 


MAKE REGIS PLOT OF DATA —' 


SET PLOT,3 

IF BQ 2 THEN !NOBRAS«l 
IPSVM-O 

IF JJ EQ 1 THEN OPBNW,5,TITLE +' COMPOSITE.PLT/Non*' 
PLOT TO,*5 “ 

IF j7 EQ 1 THEN SET VIEWPORT,.08,.48,.60,.95 
IF JJ EQ 2 THEN SET'VINWPORT, .58,.98, .60,.95 
IGRIO-O 

SCALERR - HAX(X) 

IXHAX - SCALERR 
!C-0 

PLOT,X,Y 
STOP 


■mmmmmm MAKE REGIS PLOT OF FOURIER TRANSFORM —< 
!NOERAS-l 


IF JJ EQ 1 THEN SET VIEWPORT, .08,.48,.10,.45 
IF JJ EQ 2 THEN SET^VZEWPORT,.58,.98,.10,.45 

IGRID-1 


!MTITLB-'FFT OF ' TITLE 

IXTZTLE-'FREQUENCY (MHx)' 

IYTITLB-'RELATIVE INTENSITY' 

SCCAXER - 600 ;TO CHANGE MAX FREQ CHANGE 600 TO 

; WHATB\'BR FREQ 

SCCALBR - 1200 
iXNAX- SCCALBR 

SPOT - 8CCALSR/3 
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VALUX-ABS(Z (0:92-l)) 


!C»0 

!ZGNOKX«l 

PLOT 10,8'»£Q(0:BADZE) .VALUE 
;plot,FMQ(0 :BAOZS) . VALUE 
;ir JJ EQ X TEEN STOP 

IF FMAX ME 1 THEN GOTO,FLAaOSR 
ZF JJ EQ 1 THEN BEGIN 

; aMvnaMaimnMMmiOETEBMXNE MAXIM0M8 »»m»- —» » — ««»«»—«■— a m—■ ■■■ 

READ,'USE MOUSE TO FIND MAX VALUES, BUT FIRST ENTER NUMBER OF MAXXMUM8',NMAX 
OPENN,7,TZTLE ' MAXX.OAT/LIST' 

OPENN,6,TZTLE * '"MAX.OAT/LZST' 

PRINTF,7, "S(‘^DATA VALUES OF AS)”, TITLE 
PRZNTF,6,”$('DATA VALUES OF ',AS)".TITLE 
PRINTF,7,"# INTENSITY FREQUENCY" 

PRINTF.7," 

PRINTF, 6, "^(FS .JTT'MHSR' 

CURSOR,Z1,Z2,XI,Y1 
IF XI LE 0 THEN XI - 0.0 

FOR IJ - 0,NMAX-1 DO BEGIN 
WAIT,0,5 
X2-X1 

CURSOR,Z1,Z2,XI,YX 
CORRECTO«NU 
XX-XX*X.OE6/CORRECTO 
IVAZi-MAX (VALUE (X2 :XX)) 

PRZNTF,7,'$(Z2,GX3.5,0X3.5)',ZJ,ZVAL,FRBQ(miBRX(VALUE EQ IVAL)) 
PRZNTF,6,IJ 
PRZNTF,6,IVAL 

PRZNTF, 6, FREQ (NKERB (VALUE EQ IVAL) ) 

ENDFOR 

CLOSE,6 
CLOSE,7 

ENDZF ELSE BEGIN 


/PLOTS SSNZ-LOG 
/PLOTS A LINEAR PLOT 


; M—■■■■ ■ ■■ ■■ ■■■■ ■ print MAXIMUMS to a FILE»->*»~»--«-w>----*-««-» 

READ,'USE MOUSE TO FIND MAX VALUES, BUT FIRST ENTER NUMBER OF MAXIMUMS',NMAX 
OPENN,4,TITLE + ' MAXX.OAT/LIST' 

OPENN,3,TITLE + MAX.OAT/LZST' 

PRZNTF, 4, "S('D]CTA VALUES OF ',AS)",TITLE 
PRZNTF,3,”$('DATA VALUES OF ',AS)",TITLE 
PRZNTF,4,"» INTENSITY FREQUENCY" 

PRINTF, 4," _ _" 

PRZNTF, 3, (F8 ;3T'" i HMDi 

CURSOR,ZX,Z2,XX,YX 
IF XX LS 0 THEN XX - 0.0 

FOR XJ ■ 0,NNAX-X DO BEGIN 
NAZT,0.5 
X2»X1 

CURSOR,ZX,Z2,XX,YX 
CQMUECTO-NU 
X1«X1*1.OBf/CORRECTO 
ZVAL-MAX (VALUE (X2:XX)) 

PRXNTr,4,'$(Z2,GX3.5,G13.5)',ZJ>X,:VAL,FRBQ(NHXRX(VALUS EQ ZVAL)) 

raZNTF,3,ZJ 

PRZNTF, 3, IVAL 
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PRINTF,3,FREQ(WHEW5 (VALUE EQ IVAL)) 
BNDFOR 
CLOSE,3 
CLOSE,4 

EMOELSS 

FLAGGER: TRRT-0 
ENDFOR 

AA>.03 • MIN(VALUE) 

BB».90 • SCCALER 

XFOUTS,BB,AA,STRMID(ISTIME,0,11) 

CLOSE,S 
STOP 

CLEANPLOT 

ERASE 

PLOT_TO,0 

END 
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